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Abstract 

The Mossbauer Effect (ME) is frequently used to investigate magnetically ordered systems. One 
usually assumes that the magnetic order induces a hyperfine magnetic field, Bhy per f ine , at the ME 
active site. This is the case in the ruthenates, where the temperature dependence of 
at "Ru sites tracks the temperature dependence of the ferromagnetic or antiferromagnetic order. 
However this does not happen in the rare-earth intermetallics, GdRu2 and H0RU2. Specific heat, 
magnetization, magnetic susceptibility, Mossbauer effect, and neutron diffraction have been used 
to study the nature of the magnetic order in these materials. Both materials are found to order 
ferromagnetically at 82.3 and 15.3 K, respectively. Despite the ferromagnetic order of the rare 
earth moments in both systems, there is no evidence of a correspondingly large Bh yper f ine in the 
Mossbauer spectrum at the Ru site. Instead the measured spectra consist of a narrow peak at all 
temperatures which points to the absence of magnetic order. To understand the surprising absence 
of a transferred hyperfine magnetic field, we carried out ab initio calculations which show that spin 
polarization is present only on the rare-earth site. The electron spin at the Ru sites is effectively 
unpolarized and, as a result, Bh ype rfine is very small at those sites. This occurs because the 4d Ru 
electrons form broad conduction bands rather than localized moments. These 4d conduction bands 
are polarized in the region of the Fermi energy and mediate the interaction between the localized 
rare earth moments. 

PACS numbers: 76.80. +y, 71.20.Lp, 28.20.Cz, 75.50.Cc, 75.60.Ej, 75.40.-s 
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I. INTRODUCTION 



n 

Compton and Matthias [1] found that superconductivity and ferromagnetism occur in 
Laves phase compounds containing lanthanide elements and Ru. This led to the sugges- 
tion that superconductivity and ferromagnetism could coexist in Cei-^R^Rua alloys, since 
superconductivity occurs below 6 K in CeRu2 and ferromagnetism is present in that tem- 
perature range in RRu 2 , where R is Gd or Ho|2j. The phase diagram for Cei_ x R x Ru 2 was 
investigated by Wilhelm and HillebrandQ] for R=Gd, Ho, Dy, Pr. Although no long range 
magnetic order was found in the superconducting region of the phase diagrams of these 
alloys, evidence for short-range order was found in the temperature dependence of the 155 Gd 
ME and in nuclear quadrupole resonance measurements. This occurs in Cei_ x Gd x Ru2 
for a narrow range of doping about x ~ 0.1. Fischer and Peterfg] pointed out that the 
specific heat of Cei-^Gd^R^ measured by Peter et al. pi] also showed an anomalous temper- 
ature dependence. The specific heat divided by temperature, ^, for CeRu 2 showed a sharp 
jump at the superconducting transition and a rapid fall off to zero as T — > 0, as expected. 
However the jump becomes more rounded in Cei-^Gd^R^ as x increases and £ increases as 



T 



0. The value of ^ as T — > increases as x increases from 0.05 to 0.11. This anomalous 



temperature dependence was taken as evidence of a ferromagnetic contribution to ^. The 
analysis of the data did not provide a microscopic model for the nature of this contribution 
or explain how it could coexist with superconductivity. 

Evidence for the coexistence of superconductivity and short-range ferromagnetic corre- 
lations was also found in Ceo.73Hoo.27Ru2 from the temperature dependence of the 57 Fe ME 
below ~2 K and from neutron scattering data si, [9, [l0]. Since the hyperfine coupling con- 
stant of Ru is twice that of Fejll | it was expected that there would be a magnetic field at 
the Ru nucleus of about 15 T in this material. Hyperfine magnetic fields were reported at 
the Gd site in Ce a; Gdi_ a ;Ru2|^] for x > 0.2, whose temperature dependence was consistent 
with the Curie temperature. More recently, Andoh discussed the magnetic properties of a 
number of the RRuq jl2]. 

CeRu 2 has been assigned to have the cubic Laves structure(Fd3m), although recently 
Huxley et al. 13] have shown that its symmetry is lowered due to a slight variation in the 
displacement of the Ru from their cubic Laves positions. GdRu2 and HoRu 2 are in the 
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hexagonal Laves phase structure (P63/mmc). GdRu2 has also been reported in the cubic 
Laves phase. [14] The difference in structure does not determine whether the ground state is 
superconducting or ferromagnetic since NdRu2, a ferromagnet, has the Fd3m structure. Here 
we concentrate on reconciling the ferromagnet ism found in neutron scattering, transport, 
and thermodynamic measurements with Mossbauer spectroscopy. 

Mossbauer spectroscopy is a nuclear probe of the electronic properties of systems which 
has been used to investigate magnetic order in many systems. The evidence of magnetic 
order appears in the Mossbauer spectrum as a hyperfine magnetic field induced at the 
nucleus at which the Mossbauer Effect (ME) is measured. In this way the temperature 
dependence of magnetic order has been probed by the "Ru ME in both ferromagnetic and 
antiferromagnetic ruthenates|15l. 1161 ]. The magnetic order has also been investigated using 
the 57 Fe, 193 Ir, and 155 Gd ME in the intermetallic compounds of interest here. One can 
distinguish between two cases. In the first case the B hyper f ine is found at a site on which 
there is an ordered electronic moment. This is the case of the ruthenates and a number of 
rare earth intermetallics. Using the 57 Fe ME, Wertheim and Wernick[l3] measured Bhyperfine 
values in RFe2(R=Ce, Sm, Gd, Dy, Ho, Er, Tm). Bh yP erfine is found to be ~23T in spite 
of the wide range in the size of the localized moments on the R sites. By comparison, 

^hyperfine 

investigate the magnetic structure in rare earth iron ternary intermetallics. [18[] De Graaf et 

^hyperfine 

ME from a structureless spectrum. 



the value for Bhyperfine i n ferromagnetic Fe is 32 T. The 57 Fe ME has also been used to 



al.|l9j] extracted a value for B hy p er fi n e equal to 17.5T at the Gd site in GdCu 2 using 155 Gd 



In the second case, for a non-magnetic ion in a magnetically ordered lattice, the mea- 
sured ME at this site is expected to show evidence of the magnetic order through a trans- 
ferred hyperfine magnetic field. Transferred hyperfine fields at the non-magnetic Ir site were 
measured by Atzmony et al. [^(J in RIr2(R=Pr,Nd,Sm,Gd,Tb,Dy,Ho). They found a wide 
variation from 4 T in HoIr 2 to 19 T in Gdlr 2 . Transferred hyperfine fields have also been 
measured in a rare-earth matrix doped with 1% Sn and in R 2 Sn using the 119 Sn meQ. 
For Sn doped into a rare earth, these range from -5.3 T in a Tm matrix to 23.8 T in a Gd 
matrix. The values of Bhyperfine are linear in the projection of the spin of the rare earth 
moment on its total angular momentum, (g — 1) J. In R 2 Sn, these fields range from -5.5 T 
(Er 2 Sn) to 28.9 T (Gd 2 Sn). 
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The fact that the sign of the transferred hyperfine field can change suggested that there 
is competition between different contributions which align or antialign the nuclear moment 
with the ordered electronic moment. A negative hyperfine magnetic field is antiparallel to 
the ordered electronic moment. Watson and Freeman were among the first to investigate 
the origin of the hyperfine magnetic field with large scale numerical calculations based on 
the Hartree-Fock approximation 22|, |23J . Although these calculations were limited by the 
computational capabilities then available, they demonstrated a number of the qualitative 
features of the experimental data. In particular, they showed that the ordered 4/ moment 
in Gd polarizes the electron density in the opposite direction to the ordered 4/ moment 
very close to the nucleus and in the region beyond localized 4/ bands at the edge of the 
Gd ion. They showed that the largest contribution to the hyperfine magnetic field due 
to polarization of the spin density of s electrons is the result of contributions of different 
signs from different s-shells. Their calculations also showed that the polarization of the spin 
density on neighboring sites could be opposite to that of the ordered moments leading to 
the negative hyperfine magnetic field. The sign of the transferred hyperfine field can be 
modulated by varying the lattice constant, as we describe in the section on calculations 
of the electronic properties below. In the calculations of Watson and coworkers, it was 
assumed that the 5s electrons formed the conduction band which turns out not to be the 
case in GdRu2, as we will also discuss below. In addition orbital contributions to transferred 
hyperfine magnetic field, based on the assumption that it arises from the coupling of / 
electrons on the rare earth sites with s-conduction electrons, were investigated by Dunlap 
et al. 24j Local Spin Density Approximation(LSDA) calculations were first used to calculate 
hyperfine magnetic fields in ferromagnetic 3d metals by Callaway and Wang 25, 2(|. We 
use a more recent implementation of the LSDA to calculate hyperfine fields on Ru sites and 
compare with our measured values. 

We present magnetic, transport and thermodynamic data on GdRu2, and magnetic and 
neutron diffraction data on H0RU2, showing that these are ferromagnetically ordered at low 
temperatures. However, our ME measurements of GdRu2 and H0RU2 show that Bh yper f ine at 
the Ru site is so small that, without the evidence of other experiments, one would conclude 
that there is no magnetic order. The almost complete collapse of the value of Bhy per f ine 
is an unexpected result in materials whose Curie temperatures are 82.3 K and 15.3 K, 



5 



respectively. Interestingly, t 
previously without comment 



re absence of Bhyperfine at the Ru sites in GdRu2 was noted 



27]. 



We calculate the electronic properties of these materials using a spin polarized fully 



relativistic all-electron linearized augmented plane wave method 



281 ] . ab initio bandstructure 



calculations have previously been used byother aut 
fields and electric field gradients (EFG).|29|, |30|, l3_l 



rors to determine hyperfine magnetic 



321 ] We find good agreement between 



the EFG tensor and that found from the Mossbauer spectra of CeRu2. We also find that 
the calculated B hyper f ine at the Ru sites in GdRu 2 and HoRu 2 are much smaller than those 
on Gd and Ho, consistent with the experimental results. First we present the experimental 
results. 



II. EXPERIMENTAL RESULTS 
A. Evidence of Ferromagnetism 

The transport and thermodynamic properties of GdRu 2 were investigated with a number 
of probes. The temperature, T, dependence of electrical resistance, R, and the slope dR/dT 
of a polycrystalline sample of GdRu 2 are plotted in Fig. [TJ A breaking curvature in R, 
accompanied by a sharp increase of dR/dT, occurs at ~ 86 K as temperature decreases, 
which is due to the development of an ordered state. The dc-magnetic susceptibility Xdc is 
measured from 1.9 K to 300 K at an applied magnetic field H = 50 Oe in the zero field cooled 
(ZFC) and field cooled (FC) states and the data are displayed in Fig. [2j Hysteresis in Xdc(^) 
appears at ~ 83 K. A Curie- Weiss analysis was done on the molar magnetic susceptibility 
data Xmoi given by 

_ i iv A/ 4 
Xmol -3MT-e C w)' 

where Na is the Avogadro's number, /z e fr = g(JLS)^/ J(J + 1)/j>b is the effective magnetic 
moment, g(JLS) is the Lande g-factor, k-Q is the Boltzmann's constant, and Ocw is the 
Curie- Weiss temperature. A positive 0cw ~ 85 K indicates a ferromagnetic order in GdRu 2 
and n e s ~ 7.6 fiB (Bohr magneton), which is close to the theoretical value 7.94 fiB of the 
free ion moment of Gd 3+ . This value of Ocw is consistent with the temperature at which 
hysteresis first appears. 

Arrott plots of magnetization M with respect to the internal magnetic flux density /io-^int 
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T(K) 

FIG. 1: Temperature T dependence of resistance R and the slope dR/dT of a polycrystalline 
sample of GdRu 2 . The value of the residual resistivity ratio RRR(300 K/2 K) is ~ 4.6. A curvature 
breaking which occurs in the R{T) data at ~ 80 K corresponds to the transition temperature of 
magnetic order. 

divided by M were constructed in an attempt to determine the Curie temperature Tc more 
accurately. A conventional Arrott plot consisting of M 2 vs (fi Q H int /M) isotherms, is shown 
in Fig. [3](a). In the simplest mean-field analysis of ferromagnetism, M 2 vs (fioH- mt /M) 
isotherms form a series of parallel straight lines near Tc, and the isotherm passing through 
the origin corresponds to T c . However, the M 2 -(fi H int /M) isotherms of GdRu 2 are slightly 
curved. Therefore, we applied a modified Arrott plot M 1 ^ vs (/JoHint/M) 1 ^ , where (3 ~ 0.28 
and 7 ~ 0.98 are the critical exponents based on the Arrott-Noakes equation [3^], and the 
results are plotted in Fig. 0(b) . The value of T c is ~ 82.3 K. 

Figure H](a) shows the specific heat C of GdRu2 from 2 K to 300 K. As T decreases to 
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T(K) 

FIG. 2: Temperature dependence of dc-magnetic susceptibility Xdc of GdRu2 measured at an 
applied magnetic field H = 50 Oe. Hysteresis occurs at T ~ 83 K, coinciding with the Curie 
temperature Tc determined from the modified Arrott plot [331] . vs T is shown in the inset and 
follows a Curie- Weiss behavior, which results in an effective moment /x e g- ~ 7.6/ib and a Curie- Weiss 
temperature 0cw ~ 85 K. 



~ 87 K, C starts to rise and peaks at ~ 81 K, indicating the ferromagnetic second order 
phase transition. From the C/T vs T 2 analysis between 7K and 14 K (inset of Figure H](b)), 
the values of the electronic specific coefficient 7 and the Debye temperature 0d of GdRu2 
are estimated to be 29mJ/mol-K, and 200 K, which seem reasonable in comparison with 
LaRu2's 7 ~ 41.6mJ/mol-K and Gd ~ 158.4 K. 34j After subtraction of electron and 
phonon contributions C C (T) + Ce(T) from C(T) (Fig. SJa)), the temperature dependence 
of the estimated magnetic entropy S milg is displayed in FigS^b). Estimated S majg reaches a 
saturated value of 1.5 R above Tc which is lower than the expected value of Rln8 (~ 2R). 
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24 




FIG. 3: (a) Conventional and (b) modified Arrott plots of GdRu2- The Curie temperature Tc of 
GdRu2 is ~ 82.3 K. The choice of exponents in (b) gives a closer to linear behavior in the critical 
region for intermediate fields. 



This could be due to a temperature dependent Od of GdRu 2 or our overestimate for the 
phonon contribution to the specific heat. 

A HoRu 2 sample was synthesized so that magnetic order could also be investigated using 
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y ~ 29 mJ/mol-K 2 
D ~ 200 K 



T 2 (K 2 ) 



T(K) 



300 



FIG. 4: (a) Specific heat C vs temperature T of GdRu2- The solid line is the sum of the estimated 
contributions from band electrons and phonons C C (T) + Cg(T). Inset of (b) C/T vs T 2 below 
14 K. The estimated values of electronic specific heat coefficient 7 and Debye temperature ®d are 
~ 29 mJ/mol-fT 2 and ~ 200 K, respectively, (b) Estimated magnetic entropy 5 mag of GdRu2 after 
the contributions from band electron and phonons is removed. The saturation value of S'mag is 



1.5R, less than Rln8. 
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FIG. 5: Temperature dependence of dc-magnetic susceptibility Xdc of H0RU2 measured at an 
applied magnetic field H = 100 Oe. Hysteresis occurs at T ~ 14.5 K. Inset v7 vs T fits well with 
a Curie- Weiss behavior, which results in an effective moment /i e fj ~ 10.6/iB and a Curie-Weiss 
temperature 0cw ~ 12 K. 



neutron diffraction. Gd has a large cross-section for neutron capture which makes investi- 
gations of GdRu2 with neutron diffraction impractical. Measurements of Xdc, performed on 
H0RU2 at H = 100 Oe and from 2 K to 300 K in the ZFC and FC conditions, are displayed 
in Fig. [5l Hysteresis in Xdc{T) occurs at ~ 14.5 K. The Curie- Weiss analysis (shown in the 
inset of Fig. [5]) indicates a ferromagnetic transition takes place near 0cw ~ 12 K in H0RU2 
with a fi c g ~ 10.6/is , which agrees with the theoretical value 10.6 hb of the free ion moment 
of Ho 3+ . An interesting feature in the ZFC Xdc data of H0RU2 is that there is a 20% drop of 
the Xdc values at 2 K compared to the value at 14.5 K (similarly a ~ 13% drop in the ZFC 
Xdc of GdRu 2 , Fig. [2]). 
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FIG. 6: (Color Online)High resolution neutron diffraction patterns taken (a) above and (b) below 
the magnetic phase transition at 15.3 K. The crosses indicate the observed data and the solid 
curves the calculated intensities from the structural refinements, and the difference is shown in the 
low part of the plots. The vertical lines indicate the angular positions of the diffraction lines for 
the nuclear Bragg peaks, and (b) for the magnetic (top) Bragg peaks. 
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TABLE I: Refined crystal structure parameters for H0RU2 at 25 K (first line) and 5 K(second line) 
Space group P6/3mmc(No.l94) o=5.2310(2)(A), c=8.8265(4)(A), F=209.16(2)(A 3 ). 



Atom 


X 


y 


z 


B(A 2 ) 






M(hb) 


Ho 


1/3 


2/3 


0.0664(2) 


0.36(3) 










1 /o 

1/3 


2/3 


0.0660(2) 


0.03(3) 


6.70(8) 


4.3(1) 


7.98(8) 


Rui 











0.39(3) 



















0.22(3) 








Ru 2 


0.1704(2) 


0.3407(2) 


3/4 


0.39(3) 










0.1711(2) 


0.3422(4) 


3/4 


0.22(3) 









Rp=6.16%, wRp=7M%, x 2 =0.8387 

7.12 9.02 1.428 



High resolution powder diffraction data were collected at the NCNR on the BT-1 high- 
resolution neutron powder diffractometer, using monochromatic neutrons of wavelength 
1.5403A produced by a Cu(311) monochromator. Soller collimations before and after 
the monochromator and after the sample were 15, 20, and 7 full-width-at-half-maximum 
(FWHM), respectively. Data were collected in the 9 range of 3° to 168° with a step size of 
0.05° at 25K and 5K, above and below the magnetic phase transition. Structural refinements 
were carried out using the GSAS program. [351] 

Detailed temperature dependent measurements of the magnetic order parameter were 
carried out on the BT9 triple axis spectrometer. A pyrolytic graphite (PG) (002) monochro- 
mator was employed to provide neutrons of wavelength 2.36 A, and a PG filter was used 
to suppress higher-order wavelength contaminations. Coarse collimations of 40, 48, and 40 
FWHM on BT9 were employed to maximize the intensity. A PG(002) energy analyzer was 
used in these measurements. Inelastic measurements were taken with a fixed final energy of 
14.7 meV. 

Fig. [6] shows the diffraction pattern obtained above (25K) and below (4 K) the magnetic 
transition. The diffuse background scattering in the magnetically disordered state is due to 
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paramagnetic scattering of the uncorrelated Ho moments, which decreases with increasing 
angle due to the magnetic form factor. The overall refinement fits at both temperatures are 
excellent. The structure is found to deviate slightly from the ideal phase, and the refined 
values for the crystal structure are given in Table 1. In the ground state the Ho ions exhibit 
long range ferromagnetic order, with both in-plane and c-axis components of the ordered 
moment also given in Table 1. We find an ordered moment of 7.98f8)/i£, which is in good 
agreement with the low temperature magnetization data of Andoh[l2j and somewhat smaller 
than the Curie- Weiss value obtained from the magnetization measurements. No evidence 
was found for moments on the Ru sites. The crystal and magnetic structures are shown in 
Fig. El 

The temperature dependence of the magnetic Bragg intensity, which is proportional to the 
square of the ordered magnetic moment, is shown in Fig. 8. The temperature dependence is 
typical for magnetic ordering, and the solid curve is a least-squares fit of the intensity to a 
mean field order parameter, which provides a good fit to the data with a Curie temperature 
of 15.30(4) K. 

The diffraction data indicate that the ordered moment is reduced from the free-ion value 
of 10.0/iB, which suggests that crystal field effects are important in this system. We therefore 
carried out inelastic neutron scattering measurements on BT-7 for several temperatures and 
wave vectors to search for crystal field excitations. Figure 9 shows a scan above the phase 
transition, at a temperature of 20 K. We see two clear excitations from the crystal field 
ground state, a strong one at 2.83(5) meV and a weaker excitation at 14.82(5) meV. The 
energies and intensities turn out to be quite similar to the Ho crystal field levels observed in 
the closely related Cei_ x .Ho :r Ru2 system for smaller x s[ 9], which has the cubic C-15 Laves 
structure. The rare earth site symmetry in the CeRu2 case is cubic, 43m, and the crystal 
field level scheme has been worked out in detail. For hexagonal H0RU2 the site symmetry is 
lower, 3m, but the crystal field levels look remarkably similar nevertheless. The width for the 
higher energy level is limited by the instrumental resolution, while the level at 2.83 meV has 
an observed full-width-at-half-maximum (FWHM) of 4.34(12) meV, which is much broader 
than the resolution of 1.5 meV. The width likely originates from exchange broadening. We 
note that both excitations are clearly magnetic in origin, as their intensity decreases with 
increasing wave vector, following the magnetic form factor dependence, and they decrease 
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FIG. 7: (Color Online)The hexagonal laves phase and magnetic structure of H0RU2 at 5 K. The 
ferromagnetically ordered moments on the Ho sites have a magnitude of 7.98(8)//#- 

in intensity with increasing temperature as the ground state occupancy is depleted. 

All the thermodynamic, transport, neutron diffraction measurements clearly demonstrate 
that ferromagnetic order develops in GdRu2 below 83 K and in H0RU2 below 15 K. 
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FIG. 8: (Color Online) Temperature dependence of the magnetic Bragg peak intensity, which yields 
an ordering temperature of 15.30(4) K. The intensity above the phase transition originates from the 
nuclear Bragg peak. Uncertainities are statistical in origin and represent one standard deviation. 



B. Ru Mossbauer Spectra 



Details of the source preparation of "Rh(Ru), velocity calibration, and the experimental 



setup for transmission spectroscopy are discussed in previous papers The sample of 

CeRu 2 was made with natural Ru and contained 100mg/cm 2 of Ru. Samples of GdRu 2 , 
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FIG. 9: (Color Online) Inelastic scattering observed at 20 K and a wave vector of 2.0 A -1 . The 
elastic peak has a magnetic component and nuclear incoherent scattering, and is resolution limited. 
Two crystal field excitations are observed at 2.83(5) and 14.82(5) meV. The low energy scattering 
is broad, with a width of 4.34(12) meV. The energies and intensities are quite similar to the Ho 
crystal fields observed in cubic Ceo.73Hoo.27R.u2 s[ 9] 
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FIG. 10: (Color Online) "Ru Mossbauer spectra of Ce .88Gdo.i 2 Ru2 at 2.2 K(o), 4.2 K(y), and 
17 K(*). Both the superconducting transition temperature and Curie temperature are ~ 4 
The temperature independent spectra show no evidence of a hyperfine magnetic field. 



H0RU2, and Ce a: Gdi_ a ;Ru2 were prepared with enriched 99 Ru(95%). Typical samples con- 
tained about 65mg/cm 2 of "Ru. This made it possible to measure well-resolved spectra up 
to 150K. In the analysis of the spectra the number of inequivalent Ru sites is different in 
these materials. There is one type of Ru site in the unit cell of the cubic laves structure 
(CeRu2), whereas there are two inequivalent Ru sites whose relative abundance is 1 to 3, in 
the hexagonal Laves structure (GdRu2 and H0RU2). 

CeRu2 is a superconductor with a transition temperature, T$c, of 6K. As the Gd content 
of Ce x Gdi_ I .Ru2 increases with x, Tsc falls an d goes to zero at x ~ 0.14. Previous authors 
have suggested that superconductivity and ferromagnetism could coexist close to this doping 
level. The spectrum of Ce .88Gdo.i2Ru 2 at 2.2 K, 4.2 K and 17K is shown in Fig. [TUl 
Extrapolating the data of Wilhelm and Hillebrand[3| gives values for T$c an d Tc both ~4 
K. There is no evidence in the temperature independent spectra for a hyperfine magnetic 
field, suggesting that ferromagnetism is ruled out at this doping. However, as doping with 
Gd increases, the low temperature phase is known to be ferromagnetically ordered and at 
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FIG. 11: (Color Online) 99 Ru Mossbauer spectrum of Ceo.2Gdo.sRu2 at 4.2 K. Extrapolating the 
data of Wilhelm and Hillenbrand [3], Tc f° r this material is ~60 K. This spectrum is very similar 
to that of Ceo.88Gdo.i2Ru2 except for additional asymmetry due to the increased Gd concentration. 
There is no evidence for the hyperfine magnetic field at the Ru site expected in a ferromagnetically 
ordered material. 

x = 0.2 one would expect to see an eighteen line magnetic spectrum at 4.2 K due to a large 
value of Bhyperfine field at the Ru site. Instead Fig. [TT] shows a quadrupole spectrum which is 
somewhat more asymmetric than at £=0.88 due to the increased Gd content. The absence 
of hyperfine magnetic field even in this sample suggests that the spectra of Ceo.8sGdo.i2Ru2 
do not rule out the coexistence of magnetism and superconductivity. 

The Mossbauer spectrum for GdRu2 is an almost temperature independent single peak 
between 4.2 K and 101 K except for a shift toward more positive velocities. This is shown 
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FIG. 12: (Color Online) 99 Ru Mossbauer spectrum of GdRu 2 at 4.2 K (o) and at 78 K (o) for a 
sample whose Curie temperature is 88.6 K. The 78 K spectrum(line) is scaled by 3.0 for direct 
comparison with the spectrum at 4.2 K. The only temperature dependence in the spectrum between 
4.2 K and 100 K, based also on measured spectra at 89 K and 101 K, is a continuous shift of the 
isomer shift to slightly more positive values. Again there is no evidence of a hyperfine magnetic 
field. 

in Fig. [T2] for 4.2 K and 78 K, where the 78K spectrum is scaled by 3 to compensate 
for the temperature dependence of the recoil free fraction so that a direct comparison can 
be made with the 4.2 K spectrum. The experimental FWHM of the spectrum is slightly 
broader than that of Ru powder (0.25 mm/s). The width of the spectrum is unchanged with 
temperature. This is also the case for the spectra at 89 K and 101 K (not shown). Contrary 
to the transport, magnetic, and thermodynamic data on the same sample presented in the 
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previous section, there is no evidence of a transferred hyperfme field at the Ru sites due to 
magnetic order on the Gd sites. 

The same apparent discrepancy between the Mossbauer spectrum and the experiments 
discussed above is seen in H0RU2. The spectrum of H0RU2 at 4.2 K is similar to that of 
GdRu2- It is a single peak with no evidence of splitting due to a hyperfme magnetic field 
even though the temperature is far below the sample's Curie temperature, 15. 3K. In order to 
get a more detailed picture of the properties of these materials we examined their electronic 
structure using ab initio calculations. 



III. ELECTRONIC STRUCTURE 



The lattice parameters and structure used in the calculations are those given by Compton 
and Matthias [l] for CeRu2 and by the results of the neutron diffraction measurements of 
H0RU2 at 5 K discussed. The lattice parameters and structure were measured by Compton 
and Matthias [1] for GdRu2 at room temperature. However, the lattice frequently expands 
below a ferromagnetic transition in rare earth intermetallic compounds |36j|. So we deter- 
mined lattice constants for GdRu2 at low temperatures by scaling the lattice constants given 
by Compton and Matthias and found that the total energy was minimised when their values 
were increased by ~ 1%. 



A. Density of States and Magnetic Order 



The absence of the expected large splitting in the Mossbauer spectra in GdRu 2 and HoRu 2 
points to a very small values of Bhyperfine given that the materials are ferromagnetically 
ordered. This is completely different from the ruthenates where the dependence of the 
Mossbauer spectra reflect the internal ordered field. 15|, [16( In particular, SrRuOs orders 
ferromagnetically at 163 K and the value of Bhyperfine is found to be 33 T at 4.2 K. 

We have calculated the electronic structure and the values of Bh ype rfine at the Ru and 



lanthanide sites using the Wien2k software package 281]. These calculations use the LSDA, 
which is implemented using an extension of the augmented plane wave(APW) method. Up 
to 729 inequivalent k points were used in the P63/HIHIC structure and RKMAX was set 
at 8 to ensure convergence. The energy and charge convergence variables were ec = 10~ 5 
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Rydbergs and cc = 10 _5 e. In each region a wavefunction basis set is chosen to optimize the 
calculation. The spin on each site, S, is the net electronic spin polarization in the sphere 
surrounding that site. 

The calculated contributions to the densities of states in CeRu2, GdRu2, and H0RU2 show 
that the /-bands are narrow in each case, pointing to predominantly localized states. In 
CeRu2 the calculated ground state is paramagnetic and the f-band is mostly above the Fermi 
energy, with small weight at Ep so that the / electrons also contribute to the conduction 
band. This is in agreement with the earlier ab initio calculations of Yanase 



371 ] and Higuchi 



and Hasegawa[38| and with experiment. 

The calculated contributions to the density of states of GdRu2 from the Gd 4/ band and 
from the Ad bands associated with the two inequivalent Ru sites sites are shown in Fig. [T31 
The Gd spin up /-band is ~4 eV below the Fermi level and spin down ^band is ~1.5 eV 
above. The calculated value of net spin, S, on the Gd sites is 3.46 which is very close to 
the Hund's rule result for Gd, |, and to the result of the Curie- Weiss fit to the magnetic 
susceptibility. The calculated values of S at the two inequivalent Ru sites in GdRu2 are 
very small, ~-0.06 and ~-0.09. The bottom two panels in Fig. [13] show the contributions 
to the density of states from Ru 4d electrons at the two inequivalent Ru sites. is the 
contribution for states with spin parallel to the moment of the Gd site and is positive, and 
the contribution from spin down, iV|, is negative for each of the inequivalent Ru two sites. 
At the Fermi energy the density of states is dominated by the Ru 4c? electrons, which form 
broad conduction bands. The contributions from Ru 4ci-electrons, iVj and iVj_, and iVf — N± 
are plotted in Fig. [TH in the energy range which is ~0.1 eV on either side of the Fermi 
energy, E^=0. It is seen that the rf-electrons from both Ru sites are polarized parallel to 
the Gd moment with ^ + ^ x —0.2 for each band arising from the two inequivalent Ru sites. 
The total Ru d density of states ~1.0 (eVRu) -1 at E=E F . 

On the other hand the Gd s-electrons have negligible weight at E_p, ~0.02 (eVGd) -1 and 
are unpolarized. This is also the case for Ru s-electrons. This suggests that Ru d-electrons 
mediate the interaction between localized Gd /-moments. In H0RU2, the Ru 4c? Ru electrons 
play the same role but the minority spin band is at the Fermi energy so that it is partially 
occupied. 

In H0RU2, the value of S at the Ho site is ~1.82. This value of S is slightly lower than 
2, the value consistent with the effective moment given by the fit to magnetic susceptibility, 
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as determined by Hund's Rules. The absence of moments on the Ru sites in H0RU2, S~0, 
is consistent with the neutron diffraction data. 

According to these results, the magnetic properties can be described by a model in which 
localized Gd or Ho /-moments couple to itinerant Ru <i-electrons, rather than to itinerant s 
electrons [241]. 



B. Hyperfine Magnetic Fields and Electric Field Gradients 

Bhyperfine at a given site has a number of contributions which have been derived with 
relativistic corrections by Bliigel et al.[39j and can be calculated directly using the Wien2k 



package [4dj] . These contributions are B con =B core +B va i ence , the Fermi contact term, B^p, the 
dipolar field from the on-site spin density, B or b, the field associated with the on-site orbital 
moment, and B/ Qt , the classical dipolar field from all other atoms in the system carrying 
moments. B core is the contribution due to polarization of core electrons and B m / ence is the 
contribution from the polarization of the valence or conduction band electrons. These are 
by far the largest contributions. Powder samples are used in the experiments with small 
crystallites so that each Ru nucleus sees the average of a collection of randomly orientated 
dipolar fields which presumably sum to zero. Therefore, we take Bi at to be zero in the 
current analysis. 

The calculated Bh yper f ine and local moments at the Ce and Ru sites in CeRu2 are es- 
sentially zero, consistent with the absence of magnetic order. The calculated contributions 
to Bhyperfine are shown in Table 2 for GdRu2 and in Table 3 for H0RU2. The values of 
Bhyperfine at the Ru sites are surprisingly small compared to that at the rare earth sites 
in both compounds, consistent with the Mossbauer data. By comparison, the transferred 



hyperfine field at the Ir site in Ir .oiFeo.99 is 14.3 T compared to 32 T at the Fe site 20j. The 
different contributions to Bh yP erfine on the Ru sites are all small, leading to modest values for 
Bhyperfine- The values of B hyper f ine at the Ru nuclei are sensitive to the value of the lattice 
constants. In order to determine the lattice constants for GdRu2 at low temperature, the 
electronic properties were calculated for values of the lattice constants given by Compton 
and Matthias|l] scaled by a common factor, a, from 0.995 to 1.05. Although the value of 
the net spin on the Gd site did not change for this range of lattice constants, the hyperfine 
magnetic fields at the Ru sites varied from -4.49 T at a = 1.00 to 0.1 T at a = 1.05. The 
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FIG. 13: (Color Online) Contributions to the Density of States in GdRu2- Lines are spin up and 
dashed-lines are spin down. E=0 is the Fermi level. 
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FIG. 14: (Color Online) Polarization of (i-electrons from the two inequivalent Ru sites. Black line 
- JVf — N\ from Ru site-1. Red line - JVf — N\ from Ru site-2. Blue and green lines above zero 
line are spin up and the corresponding dashed-lines below the zero line are spin down. E=0 is the 
Fermi level. 



lowest total energy for a formula unit occurs at a = 1.01 and the values for GdRu2 are 
calculated with this a. From Table 2, Bh V p e r f i ne at both Ru sites is negative, as one would 



expect from transferred hyperfine fields 
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23]. 



Examining the contributions to ^hyperfine on the rare earth site one sees that B core and 
Byaience are large and of opposite sign. This is again consistent with the original discussion 



of Watson and Freeman 
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23l |. The other contributions are negligible by comparison in 



GdRu 2 . However, in the HoRu 2 site B orb ~40 T, is a substantial fraction of B hyper f ine on the 
Ho site due to the almost complete cancellation of the B core and B va i ence . This difference in 
the value of B orb at the Gd and Ho sites is a reflection of the difference in the orbital angular 
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FIG. 15: (Color Online) Mossbauer spectrum of GdRu2 at 4.2 K (o) for a sample whose Curie 
temperature is 83.2 K. The solid line is the spectrum calculated with the electric field gradient 
components and the ^>hyperfine-, shown in Table 2 and 4. The fitting parameters are the isomer 
shifts at the two inequivalent Ru sites, ISi=0.0 mm/sec and IS2=0.12 mm/sec, and the half- width 
of the lines, T=0.07 mm/sec. 



momentum quantum number, L, on the Gd and Ho sites. Whereas Hunds' rules give L = 
on the Gd site because of the half-filled 4/ shell, they give L = 6 on the Ho site. 

The calculated components of the EFG tensor for CeRu2, GdRu2, and H0RU2 are given 
in Table 3. Since the Ce site symmetry is cubic the electric field gradient is zero. The 
components of the EFG at the Ru site, V xx =V y j / =2.80xl0 21 Vm -2 , are small compared to 
those in Ru02|27j, which also has a pure quadrupole spectrum, suggesting that the Ru 
sites have almost cubic symmetry. In GdRu2 and H0RU2, the EFG are different at the two 
inequivalent Ru sites. Again the magnitude of the components of the EFG suggests that 
there is almost cubic symmetry at these sites. 

The Mossbauer spectra for these materials can now be almost completely determined 
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FIG. 16: (Color online) 99 Ru Mossbauer spectrum of H0R112 at 4.2 K. The Curie temperature is 
15.3 K for this sample. The solid line is the spectrum calculated with the electric field gradient 
components and the ^>hyperfine, shown in Table 2 and 4. The fitting parameters are the isomer 
shifts at the two inequivalent Ru sites, ISi=0.05 mm/sec and IS2=0.18 mm/sec, and the half-width 
of the lines, T=0.08 mm/sec. 

with these calculated hyperflne magnetic fields and electric field gradients. The additional 
parameters are the half-width half-maximum, T, of the absorption lines and the isomer shifts 
at the Ru sites. These are choosen to fit the data. The GdRu2 spectrum at 4.2 K is shown 
in Figure [TS] and that of H0RU2 in Figure [TBI The line through the spectrum is determined 
by the EFG tensor and values of Bh yper f ine at the two inequivalent Ru sites calculated with 
the parameters in Table 1, Table 2, and Table 3. The non-Lorenzian shape is due to the 
calculated difference in the EFG's at the two sites. 

In Figure [T7| we show the Mossbauer spectrum of CeRu2 at the Ru site at 4.2 K. The 
full line is the spectrum determined with values of the electric field gradient tensor (EFG) 
at the Ru site in Table 3. Other parameters in these spectra are the isomer shift, IS=0.12 
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TABLE II: GdRu2: Net electron spin, S, and contributions to the hyperfine magnetic fields(T) at 
the Gd and Ru sites 





S 


Bcore 


^valence 


Borfr 


B<iip 


^hyperfine 


Gd 


3.46 


334.55 


-285.51 


-2.88 




-.29 


49.05 


Rui 


-0.06 


-4.83 


0.63 


-.08 




0.01 


-4.20 


Ru 2 


-0.09 


-3.78 


1.06 


-.0.20 


0.05 


-2.72 


TABLE III: H0RU2: Net electron spin, S, and contributions to the hyperfine magnetic fields(T) at 
the Ho and Ru sites 




S 


Bcore 


^valence 


Hor6 


H<iip 


^hyperfine 


Ho 


1.82 


179.32 


-161.31 


38.72 




-.41 


56.34 


Rui 


0.028 


-1.63 


-0.26 


-0.02 




0.01 


-1.90 


Ru 2 


-0.009 


-2.61 


0.69 


-0.75 




-0.070 


-2.74 



mm/sec, and the half-width of lines, r=0.11 mm/sec. 

The calculations discussed here have demonstrated that it is possible for transferred hy- 
perfine magnetic fields at the Ru sites to be very small in these ferromagnetic materials and 
allow us to reconcile the apparent contradiction between the results of the "Ru ME mea- 
surements and those of the transport, magnetic susceptibility, magnetization, and neutron 
diffraction experiments on GdRu 2 and H0RU2. 



TABLE IV: Diagonal elements of the Electric field gradient tensor (10 21 V/m 2 ) at different sites in 
GdRu2 and H0RU2. 



Ce 


(0.,0.,0.) 




Gd 


(-1.40,-1.40,2.79) 




Ho 


(-1.34,-1.34,2.68) 


Ru 


(2.80,2.80,-5.60) 




Rui 


(3.42,3.42,-6.83) 




Rui 


(3.85,3.85,-7.71) 








Ru 2 


(1.86,0.90,-2.76) 




Ru 2 


(1.30,0.90,-2.19) 
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FIG. 17: (Color Online) "Ru Mossbauer spectra of CeRu2 at 4.2 K(o). The line is the spectrum 
determined by the calculated EFG's, Table 3, using 0.09 mm sec -1 as the full-width half-maximum 
and isomer shift equal to 0.09 msec -1 . 

IV. DISCUSSION 

A basic assumption of using the "Ru ME is that the spectrum reflects the electronic 
environment in the material and that magnetic order is reflected in an induced hyperfine 
field. However the transferred hyperfine magnetic fields in ferromagnetic GdRu2 and H0RU2 
are so small that one would conclude that these materials do not have magnetic order. The 
calculated properties of these materials have shown how this apparent discrepancy between 
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Mossbauer, the collapse of the hyperfine magnetic field, and results of neutron diffraction, 
magnetization, transport and specific measurements on the same samples arise. Analyzing 
the calculated contributions to B hyper f ine , the net electronic spin, S, on the Gd and Ho in 
Tables 1 and 2 is roughly proportional to the magnitude of the largest contributions, B core 
and B va i ence . This suggests that B hyper f ine is ~ on the Ru sites because S is ~ on these 
sites. This arises because Ad Ru electrons form polarized conduction bands rather than 
localized moments. Consequently the "Ru ME is misleading regarding magnetic order in 
the RRU2 intermetallics. 

In contrast, moments in the magnetically ordered ruthenates are shared between Ru and 



O sites 
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42j | and the temperature dependence of magnetic order is reflected in B 



hyperfine 



determined from the Mossbauer spectra. The 2p orbitals of the six neighboring oxygen 
atoms in the RuO octahedra are strongly hybridized with Ru Ad electrons. A possibile 
explanation for the different behavior of Ru in the ruthenates and intermetallics lies in the 
difference in the electronegativities between Ru and O compared to that between Ru and 
rare earth atoms. The Pauling electronegativity of oxygen atoms is 3.44, which is larger than 
that of Ru(2.2)[43j. In GdRu2 and H0RU2, the opposite is the case. For Gd and Ho, the 
electronegativity values are 1.1 to 1.25. As a result the rare-earth atoms lose their 5d and 
two 6s electrons which take on the character of Ru d electrons. Competition for electrons 
between Ru atoms leads to the wide d-bands seen in the bandstucture. The structure can be 
pictured as a lattice of positive Gd and Ru ions, with the 4/ moments tightly-bound on the 
rare-earth sites as in the rare-earth metals 44j]. The Ad derived conduction bands mediate 
the coupling between the rare earth moments rather than forming localized moments on the 
Ru sites. 

There seem to be few other examples where hyperfine magnetic fields have not been 
induced by magnetic order. There is evidence for a collapse of B hyper f ine in hexagonal close- 
packed (hep) Fe, which is the stable phase at high pressures. Whereas there is evidence 
that this phase is antiferromagnetically ordered from Raman scattering, the six-line 57 Fe 
Mossbauer spectrum disappears at the transition from bec Fe to hep Fe with increasing 
pressure [45]. B hyper f ine was calculated by Seinle-Neumann et al. also using the wien2k 
software. They found that there was almost complete cancellation of the large core and 
valence contributions to B con as the atomic volume was reduced, simulating increasing pres- 
sure. In determining the lattice constants for GdRu2 at low temperature, we also found that 
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the sign of Bhyperfine changed as the lattice constants were varied, although not at the values 
of interest. However, the contributions to Bh yper f ine in GdRu2 and H0RU2 are less than 6 T 
due to the itinerant nature of the Ad electrons. 

The results of this investigation point to the unsuspected sensitivity of a nuclear probe, 
such as the Mossbauer Effect, to the details of the electronic structure in magnetically 
ordered materials. The ab initio calculations provide a quantitative description of how the 
collapse of the hyperfine magnetic field at the Ru site arises. In doing so, it was shown that 
the Mossbauer spectrum can be a probe of bandstructure, as well as the local electronic 
environment. 
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